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The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours 
above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise 
to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, 
neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of 
their shower signals in the water-Cherenkov detectors. In this paper we present for the first time 
an analysis based on down-going neutrinos. We describe the search procedure, the possible sources 
of background, the method to compute the exposure and the associated systematic uncertainties. 
No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. 
Assuming an E~ 2 differential energy spectrum the limit on the single flavour neutrino is E 2 dN/dE < 
1.74 x 10" 7 GeV cm" 2 s _1 sr _1 at 90% C.L. in the energy range 1 x 10 17 eV <£<lx 10 20 eV. 

PACS numbers: 95.55.Vj, 95.85.Ry, 98.70.Sa 



I. INTRODUCTION 

Neutrinos play a key role in the understanding of the 
origin of ultra high energy cosmic rays (UHECRs) . Their 
observation should open a new window to the universe 
since they can give information on regions that are oth- 
erwise hidden by large amounts of matter in the field of 
view. Moreover, neutrinos are not deviated by magnetic 
fields and would point back to their sources. 

In the EeV range, neutrinos are expected to be pro- 
duced in the same sources where UHECRs are thought to 
be accelerated, as well as during the propagation of UHE- 
CRs through the cosmic microwave background (CMB) 
radiation The latter are called cosmogenic neutri- 
nos and their presence is expected if the UHECRs above 
the spectral cut-off reported in Q contain a significant 
fraction of protons [S-Q . 

There are many current programs to search for high 
energy neutrinos with dedicated experiments Al- 
though the primary goal of the Pierre Auger Observa- 
tory Surface (SD) and Fluorescence Detectors (FD) is 
to detect UHECRs, UHE neutrinos (UHE^s) can also 
be identified and limits to the diffuse flux of UHE^s 
in the EeV range and above have been set using ear- 
lier Auger data |12hT4| | . Earth-skimming r neutrinos are 
expected to be observed through the detection of show- 
ers induced by the decay of emerging r leptons which 
are created by v r interactions in the Earth [15| . Using 
this mechanism for data collected from 1 January 2004 
until 30 April 2008, an upper limit was set: E^dN^/ 
dE v < 6tl x lO^GeVcnUV^r- 1 at 90% CL for each 
neutrino flavour [if])]. The SD of the Pierre Auger Ob- 
servatory has also been shown to be sensitive to "down- 
going" neutrinos of all flavours interacting in the atmo- 
sphere or in the mountains surrounding the SD, and in- 
ducing a shower close to the ground [lj, [lj], [3 • In this 
paper we present an analysis based on down-going neu- 
trinos and place a competitive limit on the all-flavour 



a now at University of Maryland. 
b now at Universit de Lausanne. 
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d now at Konan University, Kobe, Japan. 
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diffuse neutrino flux using data from 1 January 2004 un- 
til 31 May 2010. 

The main challenge in detecting UHE neutrinos with 
the Pierre Auger Observatory is to identify a neutrino- 
induced shower in the background of showers initiated by 
UHECRs, possibly protons or heavy nuclei [l9| and, in a 
much smaller proportion, even photons (20| . 

The identification of ^-induced showers is illustrated 
in Fig. [TJ If the incidence is nearly horizontal, "old" 
showers induced in the upper atmosphere by protons, 
nuclei or photons have a thin and flat front at ground 
level, containing only high energy muons and their ra- 
diative and decay products, concentrated within a few 
tens of nanoseconds. On the other hand, "young" show- 
ers, induced by neutrinos at a low altitude, have a thick, 
curved front with a significant electromagnetic compo- 
nent spread in time over hundreds of nanoseconds, spe- 
cially in their earlier part that traverses less atmosphere. 
In this work, to obtain an unambiguous identification of 
neutrinos, we select showers with zenith angle 9 > 75° 
and we apply criteria to ensure a deep interaction. Us- 
ing less inclined showers is in principle possible, but will 
require a better control of the various sources of back- 
ground. 

The method was tuned using data taken at the SD in 
the period from 1 January 2004 until 31 October 2007. 
A blind scan over the data collected in the remaining 
period, i.e., from 1 November 2007 until 31 May 2010 
reveals no candidates and we place a stringent limit on 
the diffuse flux of UHE neutrinos. 

For that purpose, we calculate the probability for a 
shower, produced deeply in the atmosphere, to trigger 
the SD and to be identified as a neutrino candidate. This 
probability depends on the neutrino flavour and type of 
interaction - charged current (CC) or neutral current 
(NC) - and is also a function of neutrino energy E u , inci- 
dent zenith angle 9, and atmospheric interaction depth. 
From these identification probabilities we calculate the 
exposure of the SD to deep inclined neutrino showers. 
We give an estimate of the systematic uncertainties on 
the diffuse neutrino flux limit, and discuss the implica- 
tions of our observations for models of UHE neutrino 
production. 
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FIG. 1. Pictorial representation of the different types of showers induced by protons, heavy nuclei and "down-going" (DG) as 
well as "Earth-skimming" (ES) neutrinos. The search for down-going showers initiated deep in the atmosphere is the subject 
of this work. 



II. THE PIERRE AUGER OBSERVATORY 

The Pierre Auger Observatory is a hybrid detector lo- 
cated in Malargiie, Mendoza, Argentina plj . It consists 
of an array of particle detectors |22| and a set of fluores- 
cence telescopes [23[ at four sites that provide a unique 
cross calibration capability. 

The SD is spread over a surface of ~3000 km 2 at an 
altitude of ^1400 m above sea level. This corresponds to 
an average vertical atmospheric depth above ground of 
X ground — 880 g cm~ 2 . The slant depth D is the total 
grammage traversed by a shower measured from ground 
in the direction of the incoming primary particle. In the 
flat-Earth approximation D = (X grounc i — Xi nt )/ cosO, 
where Xi nt is the interaction depth and 9 the zenith an- 
gle. For very inclined showers the curvature of the atmo- 
sphere is taken into account. 

The four fluorescence sites are located at the perime- 
ter of the surface array viewing the atmosphere above 
it [23| . In this work only data collected with the SD 
of the Pierre Auger Observatory are used to search for 
down-going neutrinos. 



A. The Surface Detector 

Since the beginning of its operation for physics analy- 
sis, in January 2004, the SD array has grown steadily and 
it has been recording an increasing amount of data. It 
consists of ^1660 detector units (water- Cherenkov sta- 
tions) regularly spaced in a triangular grid of side 1.5 
km. Each detector unit is a cylindrical polyethylene tank 
of 3.6 m diameter and 1.2 m height containing 12,000 
liters of purified water. The top surface has three pho- 
tomultiplier tubes (PMTs) in optical contact with the 
water in the tank. The PMT signals are sampled by 
flash analog digital converters (FADC) with a frequency 
of 40 MHz. Each surface detector is regularly moni- 
tored and calibrated in units of vertical equivalent muons 
(VEM) corresponding to the signal produced by a fi 



traversing the tank vertically and through its center [24j . 
The surface stations transmit information by radio links 
to the Central Data Acquisition System (CDAS) located 
in Malargiie. The PMTs, local processor, GPS receiver, 
and the radio system are powered by batteries regulated 
by solar panels. Once installed, the local stations work 
continuously without external intervention. 



B. The trigger 

A local trigger selects signals, either with a high peak 
value, or with a long duration. The second condition 
favours stations hit in the early stage of the shower devel- 
opment (moderately inclined or deeply induced showers) . 
The global trigger requires either 4 stations satisfying one 
of the conditions, or 3 stations satisfying the second one, 
in a compact configuration (see [25J for more details). 

With the complete array, the global trigger rate is 
about two events per minute, one half being actual 
shower events with median energy of 3 x 10 17 eV. 



III. SIMULATION OF NEUTRINO 
INTERACTIONS, INDUCED SHOWERS AND 
THE RESPONSE OF THE SURFACE 
DETECTOR. 

Monte Carlo simulations of neutrino-induced showers 
are used to establish identification criteria and to com- 
pute the acceptance of the SD to UHEi/s. The whole 
simulation chain is divided in three stages: 

1. High energy processes: 

• The z/-nucleon interaction is simulated with 

HERWIG (26| . 

• In the case of v T CC interactions, the r lep- 
ton propagation is simulated with a dedicated 
code and its decay (when necessary) with 

TAUOLA [27j . 
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2. The shower development in the atmosphere is pro- 
cessed by AIRES . 

3. The Surface Detector simulation is performed with 
the Offline software [Hj]. 

In the next subsections we discuss each stage in detail. 



A. Neutrino interaction 

herwig is a general-purpose event generator for high- 
energy processes, with particular emphasis on the de- 
tailed simulation of QCD parton showers. Here it is used 
to compute the fraction of the primary energy that goes 
into the hadronic vertex and to provide the secondary 
particles produced for both charged (CC) and neutral 
current (NC) interactions (see Fig. [5] for a summary of 
all the channels considered in this work). 

The energy carried by the hadronic jet is always con- 
verted into a shower which could be seen by the SD. In 
addition, the energy of the lepton produced in a CC in- 
teraction may be totally or partially visible. An electron 
is promptly converted into an electromagnetic shower. A 
r at EeV energies has a decay length of ~ 50 km and 
may decay before reaching the ground producing a sec- 
ondary shower that can be detected (so called "double 
bang" event). On the other hand, it is very unlikely that 
a high energy muon will produce a detectable shower, so 
its interaction and/or decay are not simulated. For all 
channels and neutrino flavours a set of primary v inter- 
actions is constructed from a grid of incoming neutrino 
energies, zenith angles and interaction depths. In "dou- 
ble bang" events the decay products of the r lepton are 
generated by tauola. The energies and momenta of the 
secondary particles are then injected into the program 
aires to generate the atmospheric cascade. 



B. Down-going neutrinos interacting in the 
mountains 

In addition to the interactions in the atmosphere, we 
also take into account the possibility of r neutrino inter- 
actions within the mountains around the Pierre Auger 
Observatory (mainly the Andes located to the north- 
west of the array), producing a hadronic jet and a r 
lepton. The hadronic or electromagnetic showers pro- 
duced by neutrinos of any flavour are absorbed either in 
the rock itself, or in the few ten kilometers of atmosphere 
between the mountains and the Auger array, and may be 
neglected. So only showers induced by the decay of the 
rs may be seen. In other terms, this process is exactly 
equivalent to the "Earth-skimming" mechanism, but it 
is included in this study because such showers are going 
downwards. 

The topography surrounding the SD of the Auger 
Observatory is accounted for using a digital elevation 



map |3(J. For the Auger site, the line of sight intercept- 
ing the mountains corresponds only to zenith angles very 
close to the horizon (9 > 89°). Even though the solid 
angle is much smaller than for showers with 9 > 75°, this 
mechanism is still relevant because mountains are much 
more massive. It is simulated in the same way as the 
"double bang" process, accounting in addition for energy 
loss of the t lepton in the rock [31j . 



C. Detector simulation 

To avoid excessively long computing times AIRES uses 
the standard thinning procedure [32j consisting in fol- 
lowing only some branches in the tree of interactions in 
the atmosphere. Weights are attributed to the surviving 
branches, obtaining a representative set of particles at 
any stage, especially at ground level. The first step in 
the detector response simulation is to regenerate a fair 
sample of the particles expected in each station from the 
thinned output of aires. This unthinning procedure is 
detailed in [33|]. Each particle reaching a surface detec- 
tor station is injected in the station, and the amount 
of Cherenkov light produced in water calculated with 
GEANT4 [H|. The FADC traces of the PMT signals are 
simulated using the Offline framework [Hj]. The total 
signal due to the particles entering the station, as well as 
several quantities characterizing the FADC trace which 
will be relevant for neutrino identification (see below) are 
then calculated. The local and global trigger conditions 
are applied in the same way as for real data. 



IV. INCLINED EVENT SELECTION AND 
RECONSTRUCTION 

Events occurring during periods of data acquisition in- 
stabilities :2S] are excluded. After a "trace cleaning" 
procedure removing the accidental signals (mainly atmo- 
spheric muons), the start times of the signals in the sta- 
tions are requested to be compatible with a plane shower 
front moving at speed c. If this condition is not ful- 
filled using all stations included in the global trigger, an 
iterative procedure removes stations until a satisfactory 
configuration is found with at least four stations. Other- 
wise the event is rejected. The angle between a vertical 
axis and the perpendicular direction to this plane is the 
reconstructed zenith angle 9 rec of the shower. Nearly 
horizontal showers are selected by requiring 9 rec > 75°. 
In some cases a non-inclined event, produced by detector 
fluctuations or two independent showers arriving close in 
time (less that 60 ns), may be incorrectly reconstructed 
as inclined. To remove these events we also compute the 
apparent speed of propagation of the trigger between ev- 
ery pair of stations (Va) and the average speed of the 
event ((V)), as in [16J. Genuine inclined showers have 
a "footprint" (configuration of the stations) elongated in 
the direction of arrival (left panel of Fig. The appar- 
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FIG. 2. Different types of atmospheric showers induced by neutrinos. 



ent speed of propagation of the signal, along the major 
axis of the footprint, is concentrated around the speed of 
light c. 



Under the plane front approximation, the zenith angle 
is ~ arcsin(c/ (V)). In Fig. |4]we show the distribution of 
(V) for events with 9 rec > 75° acquired between 1 Jan- 
uary 2004 and 31 October 2007. The shaded region cor- 
responds to misreconstructed or low quality events (see 
right panel of Fig. [3] for an example) . To remove these 
events we optimized a set of quality cuts using a MC sam- 
ple of 5000 regular inclined showers initiated by hadrons 
near the top of the atmosphere: (V) is required to be 
less than 0.313 m ns _1 , with a relative spread smaller 
than 0.08%. Also, the "footprint" is required to be elon- 
gated: L/W > 3, where L and W are the length and 
the width (eigenvalues of the inertia tensor, as defined 
in [HI). These cuts reject only 10% of genuine inclined 
showers. 



For events where all stations are aligned along one of 
the directions of the array, 9 rec cannot be computed and 
we rely on the average speed of the event, (V). These 
"in-line" events are of great importance since the Monte 
Carlo simulations show that low energy neutrinos (< 10 18 
eV) typically present this type of configuration in the SD. 



There is an additional requirement for events consti- 
tuted by an in-line event plus a non-aligned station (a 
non-aligned event that would become in-line by remov- 
ing just one station). This kind of spatial configuration 
is particularly prone to bad reconstruction if the non- 
aligned station was triggered by accidental muons not 
belonging to the shower front. To avoid this problem 
we also reconstruct the in-line event obtained by the re- 
moval of the non-aligned station and require it to have 
mean ground speed compatible with a zenith angle larger 
than 75°. 



V. IDENTIFICATION OF NEUTRINO 
CANDIDATES 

For this analysis, the whole data period (1 Jan 04 - 
31 May 10), was divided into two separate samples. Se- 
lected events recorded between 1 Jan 04 and 31 Oct 07 
(equivalent to ~ 1.4 years of a complete SD array work- 
ing continuosly) constitute the "training" sample, used 
to develop and optimize the neutrino identification algo- 
rithms. Data collected between 1 Nov 07 and 31 May 
10 (equivalent to ~ 2 yr of the full array), constitute 
the "search" sample. These latter events were not pro- 
cessed before the final tuning of the algorithms defining 
the neutrino identification criteria. 



A. Discrimination of neutrinos from hadronic 
showers 

Neutrinos, unlike protons and heavier nuclei, can gen- 
erate showers initiated deeply into the atmosphere. The 
main signature of these deep showers in the SD is a sig- 
nificant electromagnetic (EM) component spread in time 
over hundreds of nanoseconds, especially in the region 
on the ground at which the shower arrives earlier (see 
Fig-El- On the other hand, hadron-induced showers start 
high in the atmosphere, their electromagnetic component 
is fully absorbed and only high energy muons and their 
radiative and decay products reach the surface, concen- 
trated within a few tens of nanoseconds. 

We identify stations reached by wide EM-rich shower 
fronts via their Area-over- Peak ratio (AoP), defined as 
the ratio of the integral of the FADC trace to its peak 
value, normalized to 1 for the average signal produced 
by a single muon. In background horizontal showers the 
muons and their electromagnetic products are concen- 
trated within a short time interval, so their AoP is close 
to 1. In the first stations hit by a deep inclined shower, 
it is typically between 3 and 5 (see left panel of Fig. [5]). 

To quantify the distinctive features of hadronic and 
deeply penetrating showers induced by neutrinos at large 
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FIG. 3. Left panel: Event produced by a nearly horizontal shower (9 rec = 80°). The footprint (ellipse) is elongated along the 
reconstructed direction of arrival (arrow). Right panel: a non- inclined event with 9 rec = 79°. The major axis of the footprint 
and the reconstructed direction of arrival do not point in the same direction. Close inspection of the event suggests that stations 
3 and 5 are accidental and corrupt the reconstruction. The numbers indicate the triggering order of the stations. 
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FIG. 4. Distribution of the mean ground speed of the signal 
for events with rec > 75° acquired between 1 January 2004 
and 31 October 2007. 



zenith angle, improve the separation between the samples 
and enhance the efficiency, while keeping a simple physi- 
cal interpretation of the identification process, we choose 
a multivariate technique known as the Fisher discrimi- 
nant method [35j] . To tune it we used as a "signal" sam- 
ple the Monte Carlo simulations - exclusively composed 
of neutrino-induced showers - and as "background" the 
training sample introduced above - overwhelmingly, if 
not totally, constituted of nucleonic showers. We use real 
data to train the Fisher discrimination method, instead 
of simulations of hadronic showers, for two main reasons: 

• the composition of the primary flux is not known, 



FIG. 5. Upper panel: sketch of an inclined shower induced 
by a hadron interacting high in the atmosphere. The EM 
component is absorbed and only the muons reach the detector. 
Lower panel: deep inclined shower. Its early region has a 
significant EM component at the detector level. 



and moreover the interaction models used to sim- 
ulate hadronic showers may bias some features of 
the tail of the distributions of the observables used 
in this analysis. 

the detector simulation may not account for all pos- 
sible detector defects or fluctuations that may con- 
tribute to the background to ultra-high energy neu- 
trinos, while the real data contain all of them, in- 
cluding those which are not well known, or even 
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FIG. 6. Distributions of the AoP of the earliest station (left) and the product of the first four AoP (right) in background (real 
events in the training sample) and simulated v e CC events. There is a clear separation between both samples indicating that 
the AoP of the early stations is a good discrimination observable to be used in the Fisher method. See text for more details. 



not yet diagnosed. 

Note that, since we apply a statistical method for the dis- 
crimination, the use of real data as a background sample 
does not imply that we assume it contains no neutrinos, 
but just that, if any, they constitute a small fraction of 
the total recorded events. 

After training the Fisher method, a good discrimi- 
nation is found when using the following ten variables 
14j: the AoP of the four earliest triggered stations in 
each event, their squares, their product, and a global 
early-late asymmetry parameter of the event. We in- 
clude the square of the AoP because when the distribu- 
tion of the input variables is not gaussian the addition 
of a non-linear combination of them improves the dis- 
crimination power [36|. The product of the AoP of the 
earliest four stations in the event aims at minimizing the 
relative weight of an accidentally large AoP produced, 
for instance, by a single muon which does not belong 
to the shower front arriving at a station before or after 
the shower itself. This variable is also a very good dis- 
criminator as shown in the right panel of Fig. [5] The 
early-late asymmetry parameter is a global observable of 
the event defined as the difference between the mean AoP 
of the earliest and latest stations in the event. We have 
checked in simulations that neutrino-induced events typ- 
ically have an asymmetry parameter larger than proton 
or nucleus- induced showers Finally, the addition of 
other observables characterizing the time spread of the 
signals, such as the rise-time (between 10% and 50% of 
the integrated signal) or the fall-time (between 50% and 
90% ) , or including local observables of the stations that 
trigger last in the event, do not bring about significant 
improvements in the discrimination. 

As the shower front is broader at larger distance from 
the core for both young and old showers, the discrim- 
ination is better when splitting the samples according 



to the multiplicity N (number of selected stations). A 
Fisher discriminant was built separately for 4 < N < 6, 
7 < N < 11 and N > 12. Left panel of Fig. [7J shows the 
excellent separation achieved for events in each of the 3 
sub-samples. 

Once the Fisher discriminant J- is defined, one has to 
choose a threshold value that separates neutrino candi- 
dates from regular hadronic showers. Because the pre- 
dictions of the neutrino detection rates are very low, we 
want to keep the expected rate of background events in- 
correctly classified as neutrinos well below any detectable 
signal: in practice, we wish it to be less than one event for 
each multiplicity subsample within the expected 20-year 
lifetime of the Auger Observatory. 

The training period was used to produce a reasonable 
prediction of the background. We observe that the tail 
of the background distribution of T is consistent with 
an exponential shape. In this way, we produced a fit to 
the distribution of T for the training data in the [la, 3a] 
region, where a is the RMS of the training sample. This 
procedure is illustrated in Fig. [7J We then extrapolated 
it to find the cuts corresponding to 1 event per 1, 20 or 
100 years on the full array. The validity of the extrap- 
olation is not guaranteed, but some physical arguments 
support an exponential tail, such as the fact that showers 
produced by nuclei or protons (or even photons) have a 
distribution of X max that shows an exponential shape, 
dictated by the distribution of the primary interaction. 
The exponential model may be checked below the cut 
by comparing the actual number of events observed in 
the [3cr,4(T], [4(7, 5a], [5a, 6a] and [6a, 7a] regions, to the 
number of events predicted by extrapolating the fit done 
in the [la, 3c] region. The values are in good agreement 
as shown in Fig. [Jj For our search sample (equivalent 
to two years of full detector data) we have an estimated 
background of 0.1 events for each multiplicity class that 
add up to a total of 0.3 events with a statistical uncer- 



10 



Fisher distribution - low mult. (4 < N < 6) 



Fisher distribution - low mult. (4 < N < 6) 




2 4 
Fisher discriminant value 




-2.5 



-1.5 -1 -0.5 0.5 

Fisher discriminant value 



Fisher distribution - medium mult. (7< N< 11 




Fisher distribution - medium mult. (7< N< 11) 



-10 



5 10 

Fisher discriminant value 



in 

S10 3 




3-2-10 

Fisher discriminant value 



Fisher distribution - large mult. (12 < N) 



Fisher distribution - large mult. (1 2 < N) 



in 

=10 3 



a> lu = Training data 

U] 




-10 -5 5 10 15 20 

Fisher discriminant value 



3o 4o 
12.0 3.6 
16 6 



5o 
1.1 
1 



60 
0.3 





4 -2 

Fisher discriminant value 



FIG. 7. Left panel: distribution of the Fisher discriminant (see text for details) for events with station multiplicity 4 < N < 6 
(top), 7 < N < 11 (middle), 12 < N (bottom). Real data in the training period (1 Jan 04 - 31 Oct 07) describe the nucleonic 
background, while Monte Carlo simulated down-going neutrinos correspond to the signal. The vertical lines indicate the cut in 
the Fisher value that needs to be placed to have less than 1 event in each period of time (1 yr, 20 yr, 100 yr). Right panel: fit 
of an exponential function to the distribution of the Fisher discriminant T for the training data over the [la, 3a] interval. The 
predicted (Pred.), see text, and actual (Real) number of events are given for each of the test zones ([3a, 4a], [4a, 5a], [5a, 6a] 
and [6a, 7a] ) . 
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tainty of 30%. As we do not have at present a robust 
estimation of the background systematics we take a con- 
servative approach and do not use this value to improve 
our flux upper limit. 

As can be seen in Fig. [7] the identification cuts reject 
a small fraction of the neutrino events. Consequently, its 
choice has only a small impact on the neutrino identifi- 
cation efficiency. The neutrino-induced showers rejected 
by these cuts are those interacting far from the ground 
and similar to nucleonic-induced showers. 



VI. IDENTIFICATION EFFICIENCIES AND 
EXPOSURE 

During the data taking, the array was growing and had 
sporadic local inefficiencies. Simulations of deep inclined 
neutrino showers indicate that besides an elongated pat- 
tern on the ground they have a large longitudinal uncer- 
tainty on the core position. For these reasons we cannot 
apply (as done in the case of vertical showers [25]) a ge- 
ometrical method relying on the estimated position of 
the shower core within a triangle or hexagon of active 
stations at each time. Moreover, a shower can trigger 
the surface detector even if the core falls outside the ar- 
ray. Besides, for deep inclined showers the trigger and 
identification efficiencies depend not only on the shower 
energy and zenith angle but also on the depth of the first 
interaction. For these reasons a specific procedure was 
designed to compute the time dependent acceptance and 
the integrated exposure. 

The instantaneous status of the array is obtained from 
the trigger counting files, which respond to the modi- 
fications of the array configuration at every second. To 
avoid having to cope with an enormous number of config- 
urations, we approximate the calculation of the aperture 
by subdividing the search period in three-day intervals, 
and we select a reference array configuration to represent 
each. Once this is done we calculate the neutrino iden- 
tification efficiencies and the aperture assuming that the 
array remains unchanged during each three-day interval. 
Each reference configuration was chosen so that this ap- 
proximation, if wrong, underestimates the exposure by a 
small amount (~ 1%). 

MC generated neutrino showers produced by aires 
were randomly distributed over an extended circular area 
around the array, such that a shower with a core falling 
outside this area has no chance to trigger the array. For 
each three-day configuration the FADC traces of the ac- 
tive Cherenkov stations were simulated, the local and 
global trigger conditions were applied and the events were 
processed through the same reconstruction and identifi- 
cation algorithms as the data fSec. fV|). 

Fig. [5] shows an example of a shower that would be a 
neutrino candidate in an ideal array, placed at four ran- 
dom positions on the circular surface defined above. Two 
of the realizations are effectively recognized as neutrino 
events in the real array for that particular layout. The 




Distance [km] 

FIG. 8. An example of the result of placing the same deeply 
penetrating neutrino-induced shower at 4 different positions 
in an actual array configuration (shaded area) corresponding 
to 27 Oct 2007. The arrows indicate the azimuthal arrival di- 
rection of the shower, the dots represent the infinite ideal ar- 
ray and the circumference the extended area (see text). Solid 
symbols - either circles or squares - correspond to triggered 
stations of the simulated shower that are also on the actual 
array. Open symbols are stations that are not in the real 
array. Shower 1 is completely contained and identified as a 
neutrino. Shower 2 falls entirely outside the real array and 
it does not trigger the array. Although shower 3 triggers the 
array, it is not identified as a neutrino because the earliest 
three stations are not in the real array. Shower 4 loses some 
stations but keeps the earliest which are enough to identify 
the event as a neutrino. 



other two are either not seen, or not identified as neutri- 
nos. 

Fig. [9] shows the efficiency (fraction of events which 
pass all steps) as a function of interaction depth in the 
atmosphere for neutrinos of E v = 10 18 eV and 9 = 80°, 
in an "ideal" array without holes nor edges. There is es- 
sentially a plateau between a minimal depth (needed for 
the ^-induced shower to reach a sufficient lateral expan- 
sion) and a maximal one (such that the electromagnetic 
component is almost extinguished at ground level) . Be- 
low and above this plateau, the efficiency drops rapidly 
to zero. In other words, for a given channel and given 
values of 9 and E v , there is a slice of atmosphere above 
the array where the interactions are detected and dis- 
tinguished: the matter contained in this volume will be 
referred to as the "mass aperture" in the following. 

For each three-day period, we compute the effective 
area defined as the integral of the efficiency over core 
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FIG. 9. Example of trigger and identification efficiency as 
a function of the slant depth of the interaction above the 
ground. Notice that the Fisher discriminant neutrino selec- 
tion actually keeps most of the neutrino showers that pass the 
quality and reconstruction cuts discussed in Section iTVl 



FIG. 10. SD Exposure for our search period for down-going 
neutrino-initiated showers. The total exposure is shown as a 
full line. The exposure for individual neutrino flavours and 
interaction channels is also shown. 



position: 



A cS (E u ,9,D,t) = Je(f,Eu,9,D,t) cL4. 



(1) 



The effective mass aperture M e g(E l/ , t) is obtained by in- 
tegrating over the injection depth D and the solid angle: 



M eS (E v ,t) = 2tt // sin 9 cos 



A eS (E v ,9,D,t) d9 AD. 

To compute this integral we perform a spline interpola- 
tion on the finite three-dimensional mesh where A e e is 
determined. The total mass aperture is then obtained 
summing M c g(E u ,t) over different configurations corre- 
sponding to a certain period of time. It is defined inde- 
pendently of the ^-nucleon cross-section. 

A combined exposure can be obtained by a summation 
over the search period: 



'Mi s {E„,t) 



dt 



rn 



(3) 



The sum runs over the three neutrino flavours (with frac- 
tions a/) and the CC and NC interactions; m is the mass 
of a nucleon. Here we assume a full v T <->■ mixing, 
leading to ui l = 1 for the three flavours. 

We use the z/-nucleon cross-section given in [37j (CSS 
hereafter) to compute the reference exposure of our 
search period. It is shown in Fig. [TU] as a function of 
neutrino energy. In Table U we also give the mass aper- 
ture integrated in time for all the considered channels, 
allowing the reader to compute the exposure using dif- 
ferent cross-sections or flux models. 
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TABLE I. Effective mass aperture integrated over time for the 
search period (1 Nov 2007 to 31 May 2010) for down-going 
neutrinos of the Pierre Auger Surface Detector (in units of 
[g sr s]). 



VII. SYSTEMATIC UNCERTAINTIES 

The calculation of the mass aperture of the Auger Ob- 
servatory for neutrino showers requires the input of sev- 
eral ingredients which we have selected from amongst 
conventionally used options. Some of these choices are di- 
rectly related to the Monte Carlo simulation of the show- 
ers, i.e. generator of the first neutrino interaction, parton 
distribution function (PDF), air shower development and 
hadronic model. Others have to do with the precision 
of our knowledge of the topography of the mountains 
surrounding the Observatory, and some come from the 
limitations on the theoretical models estimating, for in- 
stance, the interaction cross-section or the r energy loss 
at high energies. By adding linearly all these contribu- 
tions, our estimate of the total systematic uncertainty on 
the exposure amounts to +22% -46%. 

In the following subsections we discuss in detail the 
dependence of the exposure on each of the above men- 
tioned choices by modifying the different ingredients one 
by one. 
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FIG. 11. Identification efficiency as function of the slant depth 
for systematic uncertainties studies. Comparison between in- 
teraction generators (herwig and pythia). The rest of the 
Monte Carlo input parameters remain the same. 
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Interaction generator 
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PYTHIA [40] 
HERWIG++ [41] 


-7% 
-7% 


PDF (gen. level) 


CTEQ06m 


MSTW [42] 


-7% 


Shower Simulator 


AIRES 


CORSIKA 6.9 [43] 


-17% 



Hadronic Model QGSJETII QGSJETI [44J +2% 

SIBYLL [45J -2% 

SIBYLL (f=0.3 EoV) -1% 

SIBYLL (E=3 EoV) -2% 

SIBYLL (9=85°) 0% 

SIBYLL (0=89°) +4% 



Thinning 



-7% 



TABLE II. Summary of the relative differences (RD) between 
the reference calculation of the exposure and the calculations 
done changing one of the ingredients of the Monte Carlo sim- 
ulations at a time. The RD were obtained for zenith angle 
#=80° and energy E=l EeV unless otherwise stated. The 
statistical uncertainty of all the relative differences is ±4%. 



A. Monte Carlo simulation of the shower 

The reference Monte Carlo neutrino sample was pro- 
duced with herwig 6.5.10 [26[ as interaction generator 
in combination with the CTEQ06m (38j parton distribu- 
tion functions, AIRES 2.8 as shower simulator (thinning 
value of 1CT 6 ) and QGSJETII.03 [H as hadronic model. 

In order to assess the influence of this particular choice 
of models on the detector aperture, independent sets of 
CC v e showers were generated at 1 EeV and 80° using 
different combinations of several interaction generators, 
PDFs, shower simulators, thinning values and hadronic 
models. We chose this particular energy and angle bin 
because it is the one that contributes the most to the 
limit. 

In Fig. [TT] we show, as an example, the detection ef- 
ficiency as a function of the slant depth when using our 
reference options (herwig) and when changing only the 
interaction generator (pythia). Since the shapes of the 
neutrino-selection efficiency curves remain similar, we 
can estimate the effect of changing the interaction gener- 
ator by computing the integral of the curves and report- 
ing the relative difference (RD) between them. The same 
procedure is applied to estimate the effect of changing 
other ingredients of the simulation. A summary of this 
RDs is given in Table [TT1 

We observe that the changes in interaction genera- 
tor, PDF, shower simulator and hadronic model brought 
about a decrease of the estimated aperture, with the 
choice of the shower simulation being the dominating ef- 
fect. On the other hand, an improvement of the relative 
thinning level causes the opposite effect. Although we 
cannot recompute the aperture for all possible alterna- 
tives of the relevant ingredients, the relative differences 
reported in Table [TT] serve as an estimate of the system- 
atic dependence of our result on the simulation options. 



For each category of potential systematic effects in Ta- 
ble UH we take the maximum observed RD as an estimate 
of the corresponding systematic uncertainty. A total sys- 
tematic uncertainty of +9% -33% on the exposure is ob- 
tained by linear addition of the maximum positive and 
negative deviations. 



B. v— nucleon cross-sections and r energy loss 

We adopted the uncertainty in the ^-nucleon cross- 
section as calculated in |37| . It translates into a ±7 % 
uncertainty in the total exposure. In any case, as men- 
tioned above, Table UJ shows the Auger mass aperture 
for down-going neutrinos which does not depend on the 
v cross-section; hence the expected neutrino event rate 
(and neutrino flux limit) can be computed as necessary 
for other models and values of the v cross-section (see 
e.g. ii). 



C. Topography 

As explained in section IIIIB1 the actual topography 
surrounding the Observatory has been taken into account 
by detailed Monte Carlo simulations which include dig- 
ital elevation maps. In principle, uncertainties due to 
different tau energy loss models should not be important 
for down-going neutrinos, but, due to the fact that the 
Pierre Auger Observatory is close to the Andes, a non 
negligible contribution to the event rate from down-going 
r neutrinos interacting in the mountains and producing 
a t lepton is expected (see Table IIII[) . The systematic 
error on the total reference exposure due to this channel 
amounts to ±6%, dominated by the uncertainties on the 
cross-section and energy loss models. 
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VIII. RESULTS AND DISCUSSION 

In this section we present the calculation of the upper 
limit to the diffuse flux of UHE^s and compare our re- 
sults to some selected model predictions and discuss the 
implications. 

A. Upper limit on the diffuse neutrino flux 

Once the multivariate algorithms and selection cuts 
defining a neutrino candidate were studied and tuned 
with the Monte Carlo simulations and the training data 
sample, we applied them to the search data sample. 
We first tested the compatibility between the shapes of 
the tails of the Fisher distributions during training and 
search periods by using an unbinned Kolmogorov hypoth- 
esis test, and found them to be in agreement with p- 
values of 0.37, 0.16 and 0.17 for the small, medium and 
large multiplicity classes, respectively. 

Wc found no candidate events in the search period (see 
Fig. IT2|). The highest test zone with events in the Fisher 
distribution of the search sample is the 6-7 sigma region. 
It has only one event and we expected 2.2 from the ex- 
ponential fit to the test sample. 

The expected number of events from a diffuse flux of 
neutrinos in a given energy range is given by: 



Search sample - low mult. (4 < N < 6) 



N, 



expected 



$(E v )S(E v )dE u , 



(4) 



where £(E U ) is our reference exposure (eq.|3]and Fig. lTOl) . 
The upper limit is derived for a differential neutrino flux 
$(E V ) = k ■ E~ 2 . Also, we assume that due to neutrino 
oscillations the diffuse flux is composed of electron, muon 
and r neutrinos in the same proportion. We expect less 
than one background event after the neutrino selection 
procedure is applied to the data sample corresponding 
to the reference exposure (see section 5). Given the un- 
certainties of this estimate, the number of background 
events will be assumed to be zero, which results in a more 
conservative upper limit. A semi-Bayesian extension [48j 
of the Feldman- Cousins approach [49[ is used to include 
the uncertainty in the exposure, giving an upper limit at 
90% CL on the integrated flux of diffuse neutrinos of: 



k < 1.74 x 10~ 7 GeV cm 



- 2 s- 1 sr- 



(5) 



The effect of including the systematics from MC, v- 
nucleon cross-sections and r energy loss is to increase 
the limit by 15%. The limit is quoted for a single neu- 
trino flavour. The relative importance of each neutrino 
flavour in the determination of the upper limit can be de- 
rived from Table IIII1 which gives the expected fractions 
of neutrinos in the selected sample according to their 
flavour and interaction channel. The largest contribu- 
tion comes from v e CC. The second largest is v T CC, due 
to double-bang interactions and the large average frac- 
tion of energy going into the shower in the decay of the 



in 

S10 3 




-2 



-10 12 

Fisher discriminant value 



Search sample - medium mult. (7< N< 11) 




2 4 

Fisher discriminant value 



Search sample - large mult. (12 < N) 




2 2 4 
Fisher discriminant value 



FIG. 12. Fisher distribution of the search sample (1 Nov 07 
- 31 May 10, 2010) for events with multiplicity 4 < N < 6 
(top), 7 < N < 11 (middle), 12 < N (bottom). No neutrino 
candidates are found. 



r lepton. Our result together with other experimental 
limits [50] is shown in Fig. [13] 

Another usual way of presenting the upper bound is 
in the less-model-dependent differential form. It assumes 
that the diffuse neutrino flux behaves as ~ 1/E 2 within 
energy bins of unity width on a natural logarithmic scale, 
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Channel 


CC 


IN L/ 


lotal 


e 


33% 


5% 


38% 




13% 


5% 


18% 


r air 


24% 


5% 


29% 


r mountains 15% 




15% 


Total 


85% 


15% 


100% 



TABLE III. Expected fractions of neutrinos in the selected 
sample according to their flavour and interaction channel (CC 
and NC). These fractions are derived assuming that electron, 
muon and r neutrinos are in the same proportion in the diffuse 
flux. 
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FIG. 13. Integrated upper limits (90% CL) from the Pierre 
Auger Observatory for a diffuse flux of down-going v in the 
period 1 Nov 2007 - 31 May 2010. For comparison up-going 
v r (1 Jan 2004 - 28 Feb 09) [3 and limits from other experi- 
ments [Hol] are also plotted. 



Single flavour neutrino differential limits (90% CL) 




v energy [eV] 



FIG. 14. Differential limits (90% CL) from the Pierre Auger 
Observatory for a diffuse flux of down-going v in the period 
1 Nov 2007 - 31 May 2010 and up-going v T (1 Jan 2004 - 
28 Feb 2009) [Hi]. For reference, two recent calculations of 
this flux are shown: "GZK-Fermi" [f| takes into account the 
Fermi-LAT constraint on the GZK cascade photons; the other 
"GZK-evolFRII" 1\ adopts a strong source evolution model 
for FR-II galaxies, assumed to be the sources of UHECRs. 



Reference 


N expected P. of obs. 


GZK-Fermi 


0.1 


0.9 


GZK-evolFRII 


0.3 


0.7 


MPR-max 


2.0 


0.1 


BBR 


0.8 


0.4 


TD-Necklaces 


0.8 


0.4 


Z-Burst 


7.8 


4 x 10" 4 



and is given by 2.AA/£{E V )E V accounting for statisti- 
cal uncertainties only and assuming no background [5l[ . 
The differential limit obtained including systematic un- 
certainties is shown in Fig. 1141 together with our previ- 
ous result on up-going v T [13J and two theoretical pre- 
dictions for cosmogenic neutrinos @, 0] ■ We observe that 
we achieve maximum sensitivity in the 0.3-10 EeV energy 
range. 



B. Model predictions 

There is a wide variety of models predicting fluxes 
of neutrinos with energies in the EeV range [lj. They 
are usually separated into three groups: cosmogenic 
neutrinos [e.g. la, [7| , neutrinos produced in accelerating 
sources [e.g. [53, [54j and neutrinos of exotic origin [e.g. 
|52| |. In all these models there are parameters with un- 
known values which can change the spectral shape and 
strength of the flux. In Table HVl we give the event rates 



TABLE IV. Expected number of events using the current ex- 
posure of down-going v measured by the Pierre Auger Obser- 
vatory for several models [r3.[7l. [52T - [5l| . The third column gives 
the probabilities of observing events given that we expect 
N. 



after folding these fluxes with our reference exposure. 

Current theoretical flux predictions for cosmogenic 
neutrinos 0, 0] seem to be out of reach of our present 
sensitivity. Concerning neutrinos produced in acceler- 
ating sources there are popular models [53l . |54| which 
predict event rates which could be detected in the next 
few years. Regarding exotic models [52(, TD-Necklaces 
will be within our sensitivity range in one or two years, 
while Z-Burst models are already strongly disfavoured. 
Note that all such 'top down' models are also tightly con- 
strained by the limits of the Pierre Auger Observatory on 
the photon fraction in UHECR [2fJ. 
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